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1. Introduction 



The idea behind the silent cosmological models, introduced and discussed by Matarrese 
et al (see Refs. [p4fl , [[| and || ), is the following. Starting from any consistent initial 
configuration for barotropic perfect fluid spacetime geometries ( Ai, g, u), there exist 
two physically different phenomena that convey information between adjacent worldlines 
within the preferred timelike reference congruence u, which is identified with the average 
4-velocity of the fluid matter source. These are sound waves and gravitational waves. 
Mathematically, in the 1 + 3 covariant representation of such models (see Refs. ||, ||, 
|T0| and [Tl| ), these dynamical interactions are represented by the spatial derivative 



source terms in the evolution equations. A careful look at the perfect fluid time 
derivative equations, obtained from the Ricci identities for u and the Bianchi identities, 
reveals that these terms are comprised of either the spatial derivatives of the fluid 
acceleration, or the spatial rotation terms (often also called the 'curls') of both the 
electric and magnetic parts of the Weyl curvature, E^ v and H^. Given the barotropic 
equation of state assumption, the fluid acceleration, and consequently the spatial 3- 
gradient of the fluid pressure, are linked to the spatial 3-gradient of the fluid's total 
energy density. Pressure 3-gradients and their spatial variations generate propagating 
sound waves within the fluid, while non-zero values of the spatial rotation of either of 
the Weyl curvature variables generically induce temporal changes in the other, usually 
interpreted as propagating gravitational waves. 

Both from a mathematical and a physical point of view, it is of interest to investigate 
the case in which these spatial derivative terms vanish, so that the resultant 1 + 3 
covariant evolution equations reduce to a coupled set of first-order ordinary differential 
equations, in terms of the comoving time derivative along the fluid flow lines. In this 
case, provided the constraint equations (involving only orthogonal spatial derivatives) 
are satisfied on an initial surface and remain satisfied, the subsequent evolution along 
individual fluid worldlines within u is decoupled, i.e., the covariant time derivatives 
decouple from the orthogonally projected spatial ones. This technically describes what 
is called the silent assumption for cosmological models; it expresses the physical idea of 
an absence of any form of waves or of gravitational induction (hence, any form of causal 
communication) propagating between the worldlines of neighbouring fluid elements. 
Thus, there is no exchange of new information (i.e., information not already coded 
in the initial data) between the fluid elements. 

To realise it, Matarrese et al assumed the fluid matter source to be irrotational 
dust, which generates a spacetime geometry of purely electric Weyl curvature, i.e., the 
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magnetic part of the Weyl curvature vanishes! : 

t^ = 0, p = 0^u M = 0, H ia , = 0. (1) 

Well-known exact solutions of the Einstein field equations fall into this category. 
Spatially inhomogeneous representatives are the dust spacetime geometries given by 
Szekeres |26[] ( discussed in a nice geometrical formulation by Goode and Wainwright 
18|] ), as well as Ellis' dust subclass of LRS class II spacetime geometries, which includes 



the Lemaitre-Tolman-Bondi model and the orthogonally spatially homogeneous (OSH) 



Kantowski-Sachs model ( see Refs. M and [pl| ) . In both examples, the Weyl curvature 



tensor is of algebraic Petrov type D (see, e.g., Ref. |1| ). A Petrov type I example of 
a silent model is provided by the OSH Bianchi Type-I dust solutions (see, e.g., Ref. 

0)- 

The question is how many other such silent solutions there are. This depends on the 
consistency of the constraint equations with the time evolution equations. A previous 
paper claimining that these equations are generically consistent with each other is, 
regrettably, wrong, as indicated in f2^j and demonstrated by Bonilla et al |§. However, 



the latter paper did not determine the set of consistent solutions. We mount a systematic 
attack on that question here. 

In the following, we first give the dynamical equations defining irrotational dust 
silent spacetimes in 1 + 3 covariant form; this is done in section || We derive a covariant 
integrability condition, and show that it is satisfied at linear level, i.e., for those silent 
models which are linearised inhomogeneous perturbations of FLRW dust universes. 
However the condition is non-trivial in the non-linear case. Thus there is a linearisation 
instability in the silent models. In section |3| the silent configurations are then formulated 
in 1+3 orthonormal frame (ONF) terminology. By the assumption of vanishing magnetic 
Weyl curvature a new constraint is generated. The consistency of this constraint with 
the remaining dynamical equations is investigated in section |. We show that a sequence 
of non-trivial differential and algebraic conditions results from repeated propagation of 
the constraints along the integral curves of u. These conditions are not identically 
satisfied in general irrotational silent spacetimes, throwing into question the further 
physical analysis of silent models, and we conjecture, in particular, that there are no 
consistent spatially inhomogeneous solutions with a Weyl curvature tensor of algebraic 
Petrov type I. Finally, the results obtained are discussed in section 



f Here and in the rest of this paper, we use the standard 1 + 3 covariant formalism, thoroughly reviewed 
in Ref. ( see also Refs. @, 0, g§, @, § and §). Sign and index conventions are employed 
according to Ref. . 
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2. 1 + 3 covariant formulation 

The silent conditions Eq. (|l|) describe irrotational and pressure-free fluid matter sources 

inducing purely electric Weyl curvature. The 1 + 3 covariant dynamical equations then 
take the form 

Time derivative equations 

= - \ Q 2 -2 a 2 (2) 

h% h v a (d p(T ) = - | - o p p o vp -E^ + l a 2 hT (3) 

h p p h\ (E pa ) = - \ fi a pv - E pv + 3 a^ p E v)p - o p a E% h pv (4) 

fi =-Qfi. (5) 

Constraint equations 

§ = h p p h\{V v a pa )-\Z p : = (df (6) 

= h" p h u a (y v E pa ) -\X P : = (C 2 )" (7) 

= V p ^a UT E T p u a := (C 3 )" (8) 

= K flu := ^^)Va(W ff )u A = [V x a}^ := (C 4 ) M , (9) 

= V := h% h° v) VptkX (V T E K a ) u x = [Vx£] F := (C 5 )^ . (10) 



Gau/? equations 



3 S^ = E^ v -\ Ba pu + a MP a% -fa 2 (11) 



3 



'R =2fi- l<d 2 + 2a 2 . (12) 

The constraints in this setting have the following implications. The spatial divergence 
of the fluid rate of shear a pv is related to the spatial 3-gradient of the fluid rate of 
expansion, Z p := h u p V u Q, Eq. (6), and, analogously, the spatial divergence of the 
electric part of the Weyl curvature E pu to the spatial 3-gradient of the fluid's total energy 
density, := /i^V^/i, Eq. (7). Both the fluid rate of shear tensor and the electric part 
of the Weyl curvature tensor share a common eigenframe, a property expressed by Eq. 
(8). This result was originally obtained by Barnes and Rowlingson [Q. Additionally, the 
constraints (9) and (10) express the condition that in silent configurations as defined by 
Eq. (HD, these tensors need to have zero spatial rotation. 

Maartens |21| has shown that the constraints for generic irrotational dust spacetime 
geometries (with non-zero H^) are consistent with each other and preserved along 
the integral curves of u, if they are satisfied at an initial instant. This result has 
been extended by van Elst [HJ to generic barotropic perfect fluid spacetimes (with 



non-zero H pv and E pv ). Equations (6) - (9) are just reduced forms of the constraints 
underlying general barotropic perfect fluid spacetime geometries. The character of Eq. 
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(10), however, is slightly different. It is a new constraint arising as a consequence of 
imposing the silent conditions, Eq. (|]). Namely, the vanishing of the magnetic part 
of the Weyl curvature, H^, results in the reduction of the general evolution equation 
for Hfj, u to a constraint, Eq. (10). This is an example of the conversion process where 
the imposed assumption that a dynamical variable remains zero throughout the time 
evolution results in a further integrability condition, additional to those generically 
obtained from the general 1 + 3 covariant dynamical equations. Consistency of Eq. (10) 
with the remaining set of equations demands that 

o = h% h\ (f ct ) = h% h\ (4r (is) 

holds throughout. Using the methods and identities developed in Ref. | 2lfl , this is 
equivalent to 

- i/i(c 4 r+|5^ , (i4) 

where 

:= h% h v > ^ TKX u x ai (V 5 £ K(7 ) - 3 h^ p h v > r]^ x u x V T [ a% E a)i 

- ' \ a^ p X a + E% Z a ] V ^ x u x . (15) 

Hence, using Eqs. (6) - (10), 

= . (16) 

The 1 + 3 covariant condition (|T6"|) may be simplified and clarified by rewriting Eq. (15) 
in terms of the irreducible parts of the spatial covariant derivatives of and E^. 
This is achieved via the complete covariant decomposition of the spatial derivative of 
any rank 2 symmetric tracefree tensor field Ap V orthogonal to u, presented in Ref. [22] : 

h v h p (\/ a A TK ) = A pvp + | hp< v h p> h K (\/ T A a ) 

-|[VxAf ( ^ pW H T , (17) 

where angle brackets enclosing indices denote the spatially projected symmetric tracefree 
part, and A pvp = A <pvp> is the divergence-free and irrotational part of the spatial 
derivative of A pu , known as the 'distortion' of A pv . By Eqs. (9) and (10), the spatial 
rotation terms in Eq. (17) vanish for A pv = a pu and A^ = E pu in irrotational silent 
models. The spatial divergence terms are determined by Eqs. (6) and (7). Then Eq. (15) 
may be reduced to a form in which all the terms are spatial divergences or distortions 
of the fluid rate of shear and the electric part of the Weyl curvature: 



B 



fi K (p e + - E k( - m a + — fx (M X + — E^ Z 



r) v)pUT u 7 



The form of Eq. (18) makes clear the separate roles of the spatial divergence terms and 
the totally symmetric tracefree parts of the derivatives. 
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Lesame et al recently reported that condition (|16| ) was identically satisfied. 
However, this is incorrect |J, and a different conclusion is obtained in what follows. On 
the other hand, considering covariant and gauge-invariant linearised perturbations of 
FLRW spacetimes in the sense of Ellis and Bruni [12| , Eq. ([16]) is identically satisfied, 
as all terms occurring in Eq. (15) are of second order. Thus, in FLRW-linearised silent 
models the constraints evolve consistently. As is seen from our subsequent analysis, Eq. 
(PS|) is not identically satisfied in the exact non-linear case. Consequently this leads to a 
linearisation instability in irrotational silent cosmological models. In other words, there 
would be consistent FLRW-linearised solutions, which do not correspond to consistent 
exact solutions. A constraint similar to Eq. (10) arises in the context of the work by 
Barnes and Rowlingson [Q, who, among the conditions (H), allowed the fluid pressure 
to be non-zero instead (in which case, of course, a configuration no longer can be called 
silent). However, in their work the integrability of that constraint was never established. 

This concludes the 1 + 3 covariant discussion of the irrotational silent models. 
To pursue a detailed constraint analysis, for which purpose we make use of algebraic 
computing facilities, it turns out to be of (relative) computational ease to describe the 
problem within the related 1 + 3 ONF framework instead. After reformulating the set 
of dynamical equations (2) - (10) in terms of 1 + 3 ONF variables, the consistency of 
Eq. (10) with the remaining set is further investigated. 



3. 1 + 3 ONF formulation 

The ONF approach employs a set of linearly independent 1-form fields { u) a } defined at 
each point of the spacetime manifold (M., g) such that the line element can be locally 
expressed as 

ds 2 = r]a b u a u; b , (19) 

where 77^ = diag [ — 1, 1, 1, 1 ] , ( y/—rj = 1 ) , i.e., a constant Minkowskian frame metric. 
The vector fields { e a } dual to the 1-forms { uj a } satisfy the relation 

(u a , e b )=6\. (20) 

We choose the timelike frame vector eo to be the 4-velocity of the fluid matter flow, 
e = u. This setup allows for a 1 + 3 split of the commutator relations as well as of the 
curvature variables and their field equations ( see Ref. || ), part of which are constituted 
by the Bianchi identities. 

The commutation functions j 0- ^ are defined by 

[e tt ,e&] = 7a&e c , (21) 

where the frame vectors act as differential operators, e a (T), on any geometrical objects 
T. The purely spatial components, 7 a /3 7 , decompose into an object a a and a symmetric 
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object n a p as follows 

1 \ r .= 2a [p 5 a l] +e MS n 5a , (22) 

where e Q/ g 7 is the totally antisymmetric 3-D permutation tensor with €123 = 1 = e 123 . 
The commutation functions with one or more indices equal to zero can be expressed in 
terms of the kinematical quantities and the quantity 

n a :=l v abcd e b -e c u d , (23) 

(where e a := V u e a ), which can be interpreted as the local angular velocity of the (to 
be chosen) spatial frame { e a }. 

Equation (8) expresses the fact that one can choose a common eigenframe for a a p 
and E a/3 , i.e., it is possible to diagonalise both tensors simultaneously. As was shown 
by Barnes and Rowlingson it follows from the on-diagonal components of both the 
./^-constraint, Eq. (9), and the evolution equation for the magnetic part, Eq. (10), that 



= n n = n 22 = n 33 . (24) 

They also showed that the eigenframe { e a } of a a p and E af3 is Fermi-transported along 
u, 

Vt a = , (25) 

a condition obtained from the vanishing off-diagonal components of the evolution 
equations of both the fluid rate of shear, Eq. (3), and the electric part, Eq. (4). 

With these conditions and bearing in mind that uj^ = 0, one has that { e a } is 
spanned by four individually hypersurface orthogonal basis fields, which implies that 
local coordinates can be found on (A4, g, u) with respect to which the metric tensor 
field g is diagonal. This result holds for Weyl curvature tensors of either algebraic Petrov 
type I or type D. 

In the following analysis, we will employ tracefree-adapted irreducible frame 
components for <r al3 and E a p, defined by (see, e.g., Ref. [0 ) 

A + :=-lA n =I(A 22 + A 33 ), A_:=f(A 22 -A 33 ). (26) 

This implies 

A 2 = \[ (A + ) 2 + (A_) 2 ] . (27) 

Given these specialisations the following set describing silent cosmological models 
according to conditions ([I]) can be derived from the general 1 + 3 ONF dynamical equa- 
tions |fPT| : 
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The commutators: 

[eo,ei] = - i(e-2(T + )ei (28) 

[e ,e 2 ] = - |(0 + ( 7 + + v / 3^-)e 2 (29) 

[eo,e 3 ] = - i(6 + <7 + - v / 3<7_)e 3 (30) 

[ei,e 2 ] = - (a 2 - n 3 i) ei + (ai + n 23 ) e 2 (31) 

[e 2 ,e 3 ] = - (a 3 - n 12 ) e 2 + (a 2 + n 31 ) e 3 (32) 

[ e 3 , ei ] = - (ai - n 23 ) e 3 + (a 3 + n i2 ) e x . (33) 

Coupled subsystem of ordinary differential evolution equations: 

e o (0) = - |e 2 -2<7 2 -|/x (34) 

e (<r + ) = - |(2 6 - a + ) a + - § (a_) 2 - E + (35) 

e (a_) = - | (6 + a + ) <r_ - £_ (36) 

e (£+) = - \ \t a + - (0 + a + ) E+ + a_ £_ (37) 

e (£_) = - | // <7_ - (0 - <7+) + cr_ (38) 

eo{fi) = - 0/i . (39) 

Note that this set implies <r_ = <^ = 0. 

Remaining decoupled system of ordinary differential evolution equations: 

eo(ai) = - | (0 + er+) ai ~ 75 ™23 (40) 

e (a 2 ) = |(-20 + cr+ + v / 3(7_)a 2 - ± (<r+ - 75 <J_) ™3i (41) 

eo(a 3 ) =|(-20 + a + - v / 3ff-)a 3 + H (7 + + 73 (7 -) ni2 ( 42 ) 

e (™2 3 ) = - I (0 + cr+) n 2 3 ~ 73 °~- °i (43) 

eofasi) = i(-29 + ( r + + v / 3(r-)n3i-|K-^ <r_) a 2 (44) 

eo(ni 2 ) = |(- 2 + a+ - v 7 ^ <j_) n 12 + \ (a + + ^ <r_) a 3 . (45) 

Tracefree part and trace of 3-Ricci curvature of spacelike 3-surfaces orthogonal to u 
(GauB equation): 

*S + = - \ [ 2 ei(ai) - e 2 (a 2 ) - e 3 (a 3 ) - 4 (n 23 ) 2 + 2 (n 31 ) 2 + 2 (n 12 ) 2 

- 3 (e 2 - 2 a 2 ) (n 3i ) + 3 (e 3 - 2 a 3 ) (n i2 ) ] 

= E + -|(0 + ( 7 + ) ( 7 + + |(^) 2 (46) 
= ^ [ e 2 (a 2 ) - e 3 (a 3 ) - 2 (n 31 ) 2 + 2 (n 12 ) 2 + 2 (d - 2 a,) (n 23 ) 

- (e 2 - 2 a 2 ) (n 3 i) - (e 3 - 2 a 3 ) (n 12 ) ] 

= - i(6-2<7+)(7_ (47) 
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= (e 2 - 2 a 2 ) (n i2 ) - (e 3 - 2 a 3 ) (n 31 ) + e 2 (a 3 ) + e 3 (a 2 ) 

+ 4n 3 ini 2 (48) 

= (e 3 - 2 a 3 ) (n 23 ) - (e 1 - 2 a x ) (n i2 ) + e 3 (a 1 ) + ei(a 3 ) 

+ 4ni 2 n 23 (49) 

= (ei - 2 ai) (n 3 i) - (e 2 - 2 a 2 ) (n 23 ) + ei(a 2 ) + e 2 (ai) 

+ 4n 23 n 31 (50) 

*R = 2 [ (2 ei - 3 ax) (ai) + (2 e 2 - 3 a 2 ) (a 2 ) + (2 e 3 - 3 a 3 ) (a 3 ) 
- (n 23 ) 2 - (n 31 ) 2 - (n 12 ) 2 ] 
= 2//- § 2 + 2 a 2 . (51) 

Equations (48) - (50) arise from the condition that the 3-Ricci curvature tensor be 
diagonal. 



The constraint equations: 



= 


(ei 


- 3ai) - y/3n 2 3 + ei(0) 




(52) 


= 


(e 2 


- 3a 2 ) (a+ + VSa-) + 3n 3i (a+ - ^ <7_) 


-2e 2 (0) 


(53) 


= 


(e 3 


- 3a 3 ) (<j + - v^o--) - 3n i2 (cr + + ^= a J) 


-2e 3 (0) 


(54) 


= 


(ei 


- 3 ai) (£?+) - V3 n 23 + \ e^) 




(55) 


= 


(e 2 


-3a 2 )(E + + VSE_)+3n 31 (E+ - ^ E_) - e 2 (/i) 


(56) 


= 


(e 3 


- 3 a 3 ) (£?+ - E_) - 3 n 12 (£+ + E„ 


) - e M 


(57) 


= 


(e 2 


- 2a 2 ) (n i2 ) + (e 3 - 2a 3 ) (ra 3 i) + e 2 (a 3 ) - 


e 3 (a 2 ) 


(58) 


= 


(63 


- 2a 3 ) (n 23 ) + (d - 2ai) (n 12 ) + e 3 (ai) - 


ei(a 3 ) 


(59) 


= 


(ei 


- 2ai) (n 3i ) + (e 2 - 2a 2 ) (n 23 ) + ei(a 2 ) - 


e 2 (ai) 


(60) 


= 


(ei 


- ai) (<7_) - V3n 23 a + 




(61) 


= 


(e 2 


- a 2 ) (cr+ - ^= cr_) + n 3 i (cr + + V^CT-) 




(62) 


= 


(63 


- a 3 ) (<r+ + 75 o"-) - "12 (o+ - >/3<x_) 




(63) 


= 


(ei 


-ai)^-) - V^n^ 




(64) 


= 


(e 2 


- a 2 ) (£+ - ^ £_) + n 3 i (£+ + VSE_) 




(65) 


= 


(e 3 


- a 3 ) (£+ + 7= £-) - m 2 (£+ -v^£-). 




(66) 



Here the following correspondences exist between the 1 + 3 ONF form and the 1 + 3 
covariant form of the constraints: Eqs. (52) - (54) correspond to Eq. (6), Eqs. (55) 
- (57) to Eq. (7), Eqs. (61) - (63) to Eq. (9), and Eqs. (64) - (66) to Eq. (10), 
respectively. Equations (58) - (60) follow from the Jacobi identities. By use of the 
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commutators (28) - (30) and re-substitution from known relations, one can show that 
the Jacobi constraints (58) - (60) as well as the GauB constraints (46) - (51) are 
preserved along u. 

4. Constraint analysis 

In 1 + 3 ONF variables, the specific silent model condition that an irrotational dust fluid 
matter source induces (purely) electric Weyl curvature of zero spatial rotation — Eq. 
(10) in 1 + 3 covariant terms — takes the form of Eqs. (64) - (66). In order for the silent 
assumption, as specified by Eq. (jl]), to lead to a consistent set of dynamical equations, 
the zero spatial rotation condition must be preserved along the integral curves of the 
preferred timelike reference congruence u. The constraints (64) - (66) are propagated 
along u by application of the commutators (28) - (30), and it is straightforward to show 
that they will be preserved, given that the conditions 

= £L ei(6) + | cr_ ea(» - 2 (ai a_ + V3n 23 o+) E + 

- 2(a 1 a + + ±n 23 o-_)E_ (67) 
= (E + - ± E„) e 2 (9) +l(a + -j § a_) e 2 (/x) 

+ 2 (a 2 - n 31 ) (<r+ + ^ O E + 

+ ^a 2 (a + -V3o-_)E_-j= 3 n 31 (a + + j= 3 a_)E_ (68) 
= (E + + E.) e 3 (6) + \ (a + + ^ a.) e 3 (/i) 
+ 2 (a 3 + n 12 ) (cr + - ^ <r_) E + 

- ^a 3 (a + + V3a-)E„-^n 12 ((T+-^(T-)E- (69) 

hold. The righthand sides of these equations are equivalent to the 1 + 3 covariant 
expression (15) derived above. In the following it is assumed that the electric part of 
the Weyl curvature is non-zero. 

4-1. Petrov type D 

If the spacetime geometry is spatially inhomogeneous and its Weyl curvature tensor is 
of algebraic Petrov type D (EL = =^> a_ = 0), then, as Barnes and Rowlingson [[I]] 
proved, it is identical to the dust models of Szekeres f2H|| , or special subcases thereof, 
and the equations are consistent. This can be seen as follows. For £L = a_ = 0, 
one obtains from Eq. (61) that n 23 = 0. Then Eq. (67) is identically satisfied, while 
Eqs. (68) and (69) yield 

= E + e 2 (6) + | a+ e 2 (/i) + 2 (a 2 - n 3 i) a + E + (70) 
= E + e 3 (6) + \ a + e 3 (/x) + 2 (a 3 + n 12 ) a + E + . (71) 
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However, in that case one can derive from the constraints that 

= e 2 (6) + (a 2 - n 3 i) a+ 
= e 3 (6) + (a 3 + n 12 ) <r+ 
= ~ e 2 (/i) + (a 2 - n 3 i) -E+ 
= | e 3 (/i) + (a 3 + ma) , 

so that Eqs. (70) and (71) are identically satisfied as well. A special subcase contained 
within the Szekeres family are Ellis' LRS class II dust models (see Refs. || and |16| ). 
Here, the further conditions = e 2 (/) = e 3 (/), = a 3 = n 12 and a 2 = n 31 apply, and 
the equations are again consistent. 

4-2. Petrov type I 

If, on the other hand, the spacetime geometry is spatially inhomogeneous and its Weyl 
curvature tensor is of algebraic Petrov type I {E- ^ => a_ ^ ), then, contrary to 



what was claimed by Lesame et al [20|, Eqs. (67) - (69) do not vanish identically, but 
constitute a new set of constraints. Of course, they are trivially satisfied, if a Petrov 
type I spacetime geometry is of OSH Bianchi Type-I ( e a (f) = 0, = a a = n a p ). 

Equations (67) - (69) can be interpreted as expressions for the spatial 3-gradients 
of the fluid rate of expansion, e a (0). Propagating them along u and re-substituting 
from known relations, one obtains algebraic expressions for the 3-gradients of the fluid's 
total energy density, e Q (/i), in the form 

e iG") = /i[ a 1 ,n 2 3,(T +1 (T- ) £ + ,£-,/j] (72) 
e 2 (/i) = 9i[ a2,nsi,(T + ,(T-,E + ,E-,fi] (73) 
e 3 (/i) = m[ a 3 ,ni2,a + ,a-,E + ,E_,fi } . (74) 

Here, /i, g\ and h\ are multivariate rational expressions of the variables indicated. Each 
individual term in the numerators therein is linear in either a a or n a p, and all terms 
in the expressions contain a factor of either a power of cr_ or a power of As an 
example, in an appendix we give the precise form of fx explicitly. 

Propagating Eqs. (72) - (74) along u and re-substituting from known relations, one 
obtains algebraic expressions for the spatial commutation functions a a in the form 

ax = n 23 f 3 [ Q,a+,a-,E + ,E^,fj l ) (75) 
«2 = n 31 g 2 [ Q,a + ,a-,E + ,E_,n ] (76) 
a 3 = ma h 2 [ Q,a + ,a^,E + ,E_,n } , (77) 



where again / 2 , g 2 and /i 2 are multivariate rational expressions of the variables indicated, 
with each individual term containing a factor of either a power of cr_ or a power of E_. 
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Finally, propagating Eqs. (75) - (77) along u and re-substituting from known 
relations, one obtains purely algebraic constraints of the form 



where / 3 , g 3 and h 3 are high-order multivariate polynomial expressions of the variables 
indicated, with each individual term containing a factor of either a power of er_ or a 
power of At this stage, in 1 + 3 covariant terms one has taken the fourth spatially 
projected covariant time derivative of the condition (10) that the electric Weyl tensor 
has vanishing spatial rotation. We emphasise that, although the conditions (78) - (80) 
are derived via tetrad methods, they reflect the covariant property that the constraints 
are not identically satisfied and do not become compatible after repeated differentiation. 
It is clear from Eqs. (78) - (80) that one can attempt to satisfy these conditions in four 
different ways (modulo a cyclic permutation of the axes of the spatial frame {e a }), 
depending on the number of non-zero n a p variables: 

(i) = n 23 = ri3i = ni2] this implies a a = 0, and it is straightforward to show that 
this case simply corresponds to the OSH dust models of Bianchi Type-I, 

(ii) = n 31 = ni2, n 23 ^ 0, which implies = a 2 = a 3 and f 3 = 0, 

(iii) rii2 = 0, n 23 7^ 7^ n 3 i, which implies 03 = and = f 3 = g 3 , 

(iv) all of n 2 3, n 31 and n 12 are non-zero, which implies = f 3 = g 3 = h 3 . 

Due to the particular structure of f 3 , g 3 and h 3 , case (iv) could be solved, e.g., by 
o~_ = -v^ E- = 0, which just reproduces the Petrov type D situation discussed above. 
The issue is to see whether other non-trivial solutions in the variables 0, a + , <j_, E + , 
E_, and \x to the highly complex algebraic conditions = f 3 = g 3 = h 3 could be found. 
Additionally this would involve satisfying Eqs. (72) - (77), and then showing that the 
time derivatives of this solution are consistent. Note that we have not concluded the set 
of time derivatives needed to prove the consistent result generically, rather we ceased 
pursuing the consistency conditions beyond Eqs. (78) - (80) because of the number of 
terms involved. Given that all of these conditions were satisfied, this would establish 
the existence of spatially inhomogeneous silent models with a Weyl curvature tensor of 
algebraic Petrov type I. Unfortunately, because of their complexity, we have been unable 
to determine if there is such a solution. However, we doubt that such a solution exists. 
Nevertheless, if the contrary was true, solutions of this type would have to arise in a 
quite different manner compared to the Szekeres and OSH Type-I solutions. 

Similarly, we have been unable to find non-trivial solutions in cases (ii) and (iii). 
So the cases (ii) - (iv) each require further investigation to establish a conclusive result. 















n 23 /s[ Q,a + ,a-,E + ,E-,fji ] 
"31 03 [ ©, o+, <?-■, E + , E_,fi] 
n 12 h 3 [ Q,a + ,a_,E + ,E_,/j,} , 



(78) 
(79) 
(80) 
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5. Conclusion 

On the basis of the analysis given in the previous section, we conclude with the following 

Conjecture: There are no spatially inhomogeneous irrotational dust silent 
models, whose Weyl curvature tensor is of algebraic Petrov type I. 

In this case, defined according to Eq. (JJ), the silent assumption for irrotational dust fluid 
matter sources would only reproduce already known classes of spatially inhomogeneous 
spacetime geometries Q; they would be very restricted. This would support the 
conjecture of Matarrese et al |24| that realistic gravitational collapse scenarios ( where 



< ) should involve non-zero magnetic Weyl curvature ( with respect to u ), ^ 0. 
This raises interesting issues about how well the Newtonian solutions in a realistic 
situation can correspond to the (more accurate) general relativistic description, because 
the magnetic part of the Weyl curvature vanishes in the Newtonian limit (see Ehlers 
and Buchert 0, confirming the view of Ellis and Dunsby fl3 ], as opposed to claims by 



Bertschinger and Hamilton p[ ). Furthermore, there is a linearisation instability within 
a purely general relativistic approach to irrotational silent cosmological models, since, 
as we showed in section |l the FLRW-linearised silent models are consistent. 

We note that analysis of the kind presented here can be conducted either via 
covariant or tetrad methods. In some cases, covariant methods succeed in the analysis 
of consistency (e.g. Ref. plf), in some cases a tetrad analysis seems to be required (e.g. 
Ref. pcifl), and there are cases where either formalism can be used (e.g. Refs. || and 
[ E7J ) . We emphasise that whichever formalism is used, this does not alter the underlying 
covariant nature of the problem. 

Two final remarks should be made. First, the silent criterion for barotropic perfect 



fluids, as introduced by Matarrese et al [24, m , and as applied in this paper, demanded 



that in mathematical terms the evolution equations within the set of 1 + 3 covariant 
dynamical equations reduce to a coupled set of ordinary differential equations. However, 
it should be pointed out that a coupled set of ordinary differential equations describing 
evolutionary behaviour of relativistic cosmological models can also be obtained from 
the 1 + 3 ONF dynamical equations. The best-known example of this kind are the OSH 



perfect fluid models as discussed by, e.g., Ellis and MacCallum [0 (see also Ref. |28| ). 
Here the reduction of the evolution equations to a set of ordinary differential equations 
is achieved by choosing a 1 + 3 ONF { e a } that is invariant under the motions induced 
by the G3 isometry group of spacelike translations. Hence, according to the ordinary 
differential equations criterion, the OSH perfect fluid models can be called 1 + 3 ONF 
silent, but, in general, not 1 + 3 covariantly silent. Most of the OSH perfect fluid models 
have non-zero electric and magnetic parts of the Weyl curvature, E^ v 7^ 7^ H^ u , and, 
more importantly, non-zero spatial rotation terms thereof, ^ ^ J^ v . However, 
they do not represent wave-like solutions that convey information from one worldline to 
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another. It is not clear at present how to characterise this behaviour in a 1 + 3 covariant 
manner. 

Second, the demand that H^ u = and = goes beyond what is needed to 
establish a silent model in the 1 + 3 covariant approach to perfect fluid spacetime 
geometries, for all that is required (in the spatially inhomogeneous case) is that the 
pressure p and the spatial rotations, of E^ u and of H^, should vanish. The case 
([]]) was considered for simplicity (and even here, a full solution is not yet available). 
At present we do not know, if there are silent models with E^ v ^ ^ H^ v , but 
= I = J^, or if non-zero vorticity will appreciably broaden the class of allowed 
silent models. Preliminary investigations in the latter case seem to indicate that rotating 
silent models are tied to severe restrictions and would therefore be quite rare as well. 
However, both questions need to be clarified. 
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Note added in proof: It has been drawn to our attention that the conjecture in Sec. 5 has 
been independently arrived at by CF Sopuerta [25|, who gives an interesting analysis of 
silent models from a different viewpoint. 



Appendix 

For illustrative purposes, in this appendix we explicitly give the exact form of Eq. (72). 
We have that 

ei(/i) = A[ ai,n w ,<T+,<r-,E+,E-,f* ] 
= 4 [ a x A x + ^ n 23 B x ] / C x , 

where 

A x := 6 a + o_ E + E_ - (<r_) 2 fj, E + + ((x_) 2 {E + f 

-(a_) 2 (£_) 2 + 2£+ (£_) 2 
B x := 3 (ct+) 2 fi + 24 (a + f E + E^ - 6 a+ <x_ fiE+ + 6 a + <r_ (E + f 

- 6 (7+ a- (E-) 2 — (cr_) 2 /i E- + 4 (cr-) 2 E + EL. 

+ 6 (E+) 2 E_ + 2 {E^f 
C x : = 3 a + a_ E_ - (ct_) 2 // + (<x_) 2 E + + 2 (E_) 2 . 
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The right hand sides of the remaining conditions at this and all subsequent levels of 
differentiation can be obtained from a transformation rule related to a cyclic permutation 
of the axes of the spatial frame { e a }. This rule is given by making the substitutions 

1 (23) -> 2 (31) -> 3 (12) -> 1 (23) 

A_ -> \ (V3A+ - A.) -> - | (V3 + A_) -> A_ . 
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